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Abstract 

In this paper has been considered probability-one global convergence of NFPH 
(Newton-Fixed Point Homotopy) algorithm for system of nonlinear equations and has 
been proposed a probability-one homotopy algorithm to solve a regularized smoothing 
equation for NCP with generalized monotonicity. Our results provide a theoretical 
basis to develop a new computational method for nonlinear equation systems and 
complementarity problems. Some preliminary numerical experiments shows that our 
NFPH method is useful and promissing for difficult nonlinear problems. 
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1 Introduction 

NCP (Nonlinear Complementarity Problem) is to find x G R n such that 

Xi > 0, fi(x)>0, Xifi(x) = 0, i = l,...,n. (1.1) 

Zhao and Lip] studied several properties of a homotopy solution path associated with non- 
linear quasi-monotone complementarity problems. They established a sufficient condition to 
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assure the existence and boundedness of this homotopy solution path. Their results provide 
a theoretical basis to develop a new computational method for quasi-monotone complemen- 
tarity problems. In [2], Billups and Watson considered probability-one global convergence 
of an interior FPH (Fixed Point Homotopy) algorithm for bounded MCP(Mixed Compe- 
mentarity Problem). Their idea is to reformulate the MCP as a system of equations using 
FB(Fischer-Burmeister)- NCP function and then solve smooth approximations of this sys- 
tem with a homotopy method. Billups [3] has considered probability-one global convergence 
of FPH algorithm using smoothing of FB-NCP function for MCP satisfying a coercivity 
or generalized monotonicity and strict feasibility. Watson [3] has considered probability-one 
global convergence of FPH algorithm for monotone complementarity problem. His method 
involved reformulating the NCP as a system of smooth equations and applying a homo- 
topy method to solve this system. In the context of Newton-based methods, such smooth 
reformulations of complementarity problems are inferior to nonsmooth reformulations due 
to slow local convergence for degenerate solutions. In contrast, nonsmooth reformulations 
allow much faster (superlinear or quadratic) convergence to degenerate solutions. In [5], 
Hotta and Yoshise considered global convergence of a non-interior homotopy algorithm using 
CHKS(Chen-Harker-Kanzow-Smale)-smoothing function in case that map / is Pcr m apping 
and has an interior feasible point, or / is monotone and has an interior feasible point. In 
this paper, we propose NFPH (Newton-Fixed Point Homotopy) algorithm for nonlinear sys- 
tem and consider its probability-one global convergence, and extend it to solve a regularized 
CHKS-smoothing of nonsmooth reformulation for NCP with generalized monotonicity. 



2 Homotopy method for nonlinear system 

2.1 Homotopy map and global convergence 

Our probability-one homotopy algorithm is based on the following Theorem. 

Theorem 2.1. ([2] and [3]) Let F : R n — > R n be a C 2 -function and suppose there exists 
C 2 -map p: R m x [0, l)xi?M R n such that 

(i) the n x (m + l + n) Jacobian matrix Dp(a, A, x) has rank n on the set 

p-\0) = {( a , A, x) £ R m x [0, 1) x R n \p(a, A, x) = 0} 

(ii) for any fixed a £ R m and A = , the equation p a (X,x) = p(a, A, x) — has a unique 
solution x a £ R n , 

(Hi) p(a, l,x) — F(x) for any fixed a £ R m , 
(iv) p _1 (0) is bounded for any fixed a £ R m . 

Then for almost all a £ R m (in the sense of Lebesgue measure) there exists a zero curve 
la of p a , along which the Jacobian matrix Dp a has rank n , emanating from (0,x a ) and 
reaching a zero x of F at A = 1 . Moreover, 7 a does not intersect itself and is disjoint from 
any other zeros of p a . 
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The expression "reaching a zero" means that there exists a sequence of points {(Xk,x k )} 
in 7 a , accumulating at (l,x). The popular homotopy used often in practice is the FPH 
defined by 

p(a, X, x) = XF(x) + (1 — X)(x — a). 
In this paper, we consider NFPH defined by 

p{a, X, x) = XF{x) + (1 - X)G{x, a), (2.1) 

where G(x, a) = F(x) — F{a) + A(x — a) and A is a symmetric and positive definite matrix. 
In what follows all of consideration will be made under the following assumption: 

Assumption 1. F is a C 2 -mapping and there exists a symmetric and positive definite 
matrix A such that F'(x) + A is nonsingular for every x G R n . 



Remark 2.1. Under the assumption 1, the map p(a, X,x) defined by (2.1) satisfies conditions 



(i)~(in) of Theorem 2.1 



The following Theorem provides a sufficient condition for condition (iv) of Theorem 2.1 
to be satisfied. 

Theorem 2.2. Suppose that F : R n — > R n satisfies the assumption 1 and there exists x G R n 
and M > such that 

\\x\\ A i < M, (2.2) 
where A is such as in the assumption 1 and H^IL* — Vx T Ax. If it holds 

(x - x) T F(x) > (2.3) 

for every x G R n such that \\x — x\\ i > 2M , then each zero curve 7 a of p a (X,x) = p(a, X,x) 
defined by (2.1) is bounded for every a G R n such that 

a + A^Fia) G B = {x G R n \ \\x\\ A $ < M} 

and there exists a zero curve 7 a , emanating from (0, a) and reaching a zero x of F at X = 1 
for almost all a G R n with a + A~ l F(a) G B. In particular, if F'(x) is nonsingular, 7 a has 
finite arc length. 

Proof. Let \\x\\ i > 3M . Then we have 

lb -all i > IblLi - lla'ILi >3M-M = 2M 

II 11^45 — II ll^ 2 II ll^4 2 — 



for every a' G B , and since \\x — a'|L^ < ||^IL 2 + ll -'!!^ < %M by (2.2) , we have 
(x — x) T A(x — a) — (x — a' + a' — x) T A(x — a) = 

II 'Il2 /~ r\T \t i\ ^ ii / 1 1 2 ii ~ /ii ii /ii 

= \\x — all i — [x — a) A[x — a > \\x — a i — \\x — a \\ i \\x — a \\ i = 

II 11 A3 v ' y / — II 11^2 II 11^2 11 11 A3 



\x 



\x — a 



i 

A3 



X 



> 



using generalized Schwartz inequality. Consequently, we have 

(x - x) T A(x - a') > 



(2.4) 



x \\a — • Thus, if a' = a + A G -B , we have 



for every a' E B and x G -R n with 

(x - x) T (AF(x) + (1 - A)G(x, a)) > for x G i? n with ||x|| i > 3M, i.e. 



(x - x) T p a (X,x) > 



(2.5) 



for any A G [0, 1] by (2.3) and (2.4) because ||x — x\\ i > 2M . The inequality (2.5) implies 
that p a (X,x) 7^ for AG [0, 1] and x G {x G i? n |||x||^i > 3M} . Therefore, zero curve 7 a of 
p a is contained in set [0, 1] x {x G i? n |||x||^i < 3M} and is bounded, which proves the first 
proposition of Theorem together with Theorem 2.1 If F'(x) is nonsingular, , which implies 
that 7 a has finite arc length. □ 

Definition 2.1. A function F : R n — > R n is said to be pseudo-monotone at x if (x — 
x) T F(x) > for every x implies that (x — x) T F(x) > for every x . 



Corollary 2.1. Suppose that the assumption 1 is satisfied and solution set of F(x) = is 
bounded, i.e. i < M for every solution x of F(x) = . If F(x) is pseudo-monotone at 
each solution x , then there exists a zero curve 7 a of p a , emanating from (0, a) and reaching 
a zero x of p a at X = 1 for almost every a G R n with a + A _1 F(a) G B , where A is such as 
in the assumption 1. 



have (x — x^ 7 



and for every x G R n with \\x — x\\ i > 2M 



Proof. For every solution x of F(x) 

F(x) = . Hence, by the pseudo-monotonicity of F(x) at x , we have 



we 



x-x) T F(x) > 



(2.6) 



for every x G R n with \\x - 
the proof of the corollary. 



x ||^i > 2M , which implies (2.3). Thus, Theorem 



2.2 



completes 

□ 



Remark 2.2. Suppose that x is a solution of F(x) = 0. If F is monotone in x , F is 
pseudo-monotone at x . If F is pseudo-monotone at x, F satisfies (2.3). 



2.2 Tracking the zero curve 

The zero curve can be tracked by the procedure considered in section 3.2 of [2]. As discussed 
in section 2.1, the zero curve can be parameterized by arc length : Let (A(s),x(s)) be the 
point on 7 a of arc length s away from (0, a). Tracking the zero curve involves generating a 
sequence of points {y k } G R n+1 , with y° = (0, a), that lie approximately on the curve in order 
of increasing arc length. That is, y k = (X(si~),x(sj t )) , where {s^} is some increasing sequence 
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of arc lengths. The subroutine STEPNX from HOMPACK90 [6] is used to handle the curve 
tracking. At each iteration, STEPNX uses a predictor-corrector algorithm to generate the 
next point on the curve. The prediction phase requires the corresponding unit tangent vector 
to the curve, (y') k = (A'(s&), x'(sk)) for each iterate y k . This is accomplished by finding an 
element rj of the null space of V p a (y k ) and setting (y') k = T)/\\t)\\ or (y') k = —T]/\\r]\\, where 
the sign is chosen so that (y') k makes an acute angle with (y') k ~ l , for k > 0. On the first 
iterate, the sign is chosen so that the first component (corresponding to A) of (y')° is positive. 
At each iteration after the first, STEPNX approximates the zero curve with a Hermite cubic 
polynomial c k (s) , which is constructed using the last two points y k_1 and y k , along with the 
associated unit tangent vectors {y') k ~ l and (y') k . A step of length h along this cubic yields 
the predicted point w '° = c k (s + h) . The first iteration uses a linear predictor instead, 
which is constructed using the starting point y° and its associated unit tangent vector. 

Once the predicted point is calculated, a normal flow corrector algorithm [6] is used 
to return to the zero curve. Starting with the initial point w k, ° , the corrector iterates 
w , j — 1; 2, ... , are calculated via the formula w k ^ +l = to fej + z k, \j = 0, 1, ... , where the 
step is the unique minimum-norm . 

The corrector algorithm terminates when one of the following conditions is satisfied: the 
normalized correction step z k '^/(l + Hu^H) is sufficiently small, some maximum number of 
iterations (6 in our experiments) is exceeded, or a rank-deficient Jacobian matrix is encoun- 
tered in the Newton equation. In the first case, set y k+1 = w J , calculate an optimal step 
size h for the next iteration, and proceed to the next prediction step. In the second case, 
discard the point and return to the prediction phase, using a smaller step size if possible; 
otherwise, terminate curve tracking with an error return. In the third case, terminate the 
curve tracking, since rankV p a < n should theoretically not happen and indicates serious 
difficulty The step size h is also never reduced beyond relative machine precision. 

Finally, to emphasize the robustness and effectiveness of the proposed homotopy method, 
consider the one-dimensional equation f(x) = 0, where 

f(x) = arctan(100x)/7r + sin(5x/(a; 2 + 0.2))/2 + O.lx. 

solution to the Newton equation . This function has a unique root at x — 0. For algorithms 
that rely on descent of a merit function, this root is difficult to find because the global 
minimum of the merit function 9{x) = is in a very narrow valley (Fig.l). Nevertheless, 
the probability-one homotopy algorithm of [2] easily found the root, tracking the homotopy 
zero curve in 32 steps from a starting point x° = 0.5 . 



Our Newton-fixed point homotopy method with A = I in G(x, a) of (2.1 ) found the root 
in 2 iterations(6 steps) and 3 iterations(13 steps) from the same starting point with step size 
h = 0.63 and h — 1, respectively. As a comparison, PATH version 4.0 [7] was used from the 
same starting point. After 449 iterations, PATH terminated at x =0.24233, corresponding to 
a local minimum of 9. This function 8(x), while artificial, is representative of merit functions 
encountered in applications such as protein folding, analog circuit simulation, and aircraft 
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configuration design. This simple numerical experiment shows that our NFPH algorithm is 
a promising method. 




-2 -1.5 -1 



Figure 1: Graph of merit function 9{x) 
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3 Homotopy algorithm for NCP 

Definition 3.1. If a function if : R 2 — > R is such that <p(a, b) = <^ a > 0, b > 0, ab = , 

then f is called NCP function. 



The function ip(a, b) = a + b— y^(a — b) 2 is a popular NCP function called min-function. 
Then NCP (1.1) is equivalent to solve a system of nonlinear equations 

<p(xi,fi(x)) = 0, i = l,...,n 

and the problem (1.1) is reduced to solve the following nonlinear system . 

f(x) - y 



F(x,y) 



(3.1) 
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where y) = (f(xi, fi(x)), ...,<p(x n , f n (x))) T . The min-function ip(a,b) is nondifferen- 
tiable in case of a = b and this function is approximated by the following smooth function: 



ip^a, b) = a + b — yjfl — b) 2 + 4/x 2 



Let's approximate (3.1) by the following smooth system using (3.2): 

F»{x,y) = 



f(x) - y + fix 
^fj.(x,y) + \iy 



(3.2) 



(3.3) 



where $n(x,y) = (<f fJ ,(xi,f 1 (x)),...,(p IJi (x n ,f n (x))) T . Then F°(x,y) = F(x,y). The smooth 



approximation often used for (3.1) is 



f{x) - y 



and (3.3) is a regularized smoothing for (3.1). 



Lemma 3.1. For every \i > 0, x and y, we have 

\\$^x,y) - $(x,y)\\ < 2ny/n. 
Proof. When £ > 0, r\ > 0, we have a/£ + r\ > and 

(Vl+V- V / ^) 2 = e + ^-2 v ^T^v / ^ + ^<2£ + r/-2 v ^ v ^ = r / . 
Thus, letting ^ = (scj — yi) 2 , r]i = (2/i) 2 , i — 1, n , it holds 

||$„0r,j/) - $ (x,i/)|| 2 = ^[ v /(x i -y i ) 2 + V- v 7 ^-!/,) 2 ] 2 < W 

i=i 



which gives us the inequality (3.4). 



Let's construct a homotopy to solve ( 3.3[ ) as follows. 

//(A) = 0(1 - A), AG [0,1], /3>0, 
p B (A,*) = Ai^C*) + (1 - A)G^(z,a), 

where 

z=^* y j, G^ x) (z, a) = F»M(z) - i^ (A) (a) + A(z - a) 
Let's denote zero set of p a for fixed a by 

p-\0) = {(\,z)\ Pa (\,z) = 0, AG [0,1)} 



(3.4) 



□ 

(3.5) 
(3.6) 
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Theorem 3.1. Sup pose that the assumption 1 is satisfied for some matrix A = cl and 
F^ x '(z) defined by (3.3), where I is an identity matrix and c is a positive constant. If 



a i < M, 

1 "A? ' 



1511 i < M, M > 2Jch and 

I 11^2 — V 



z - zyF(z) > 



(3.7) 



for every z G R with \\z — z\\ i > 2M , then there exists a zero curve 7 a of p a defined by 
(3.6), emanating from (0, a) and reaching a zero z of F at A = 1 for almost all a G R 2n with 
a + A- l F^ x \a) G B = {z G R 2n \\\z\\ A i < M}. 



Proof. First, let's prove that 



[z - z" T 



F^ x \z) > 



for every z G R 2n with \\z — z\\ i > 2M. Let /i denote /i(A) for simplicity. We have 
(z - z) T F»(z) = (z- z) T F(z) + (z- z) T (F»(z) - F(z)) = 



(3. 



= {z- z) 1 F(z) + (x- x) 1 fix + (y- y) 1 ($„(z) - $(z)) + (y - y) 1 fiy 

= {z- z) T F(z) + /i(z - zfz + (y- yf^^z) - $(*)) 
= {z- ~z) T F(z) +n\\z- ~z\\ 2 + »{z - ~z) T ~z + (y- y) T {^{z) - $(*)) 
>{z- z) T F{z) + fi\\z - ~z\\ 2 - fx\\z - z\\ \\z\\ - \\y - y\\ \\%{z) - $(z) \\ 
>(z- z fF(z) + fz\\z - z\\ 2 - fi\\z - z\\ \\z\\ -\\z- z\\ \\^(z) - 
which gives us 

(z - i) T F^(z) >{z- z) T F(z) +fi\\z- z\\ 2 - n\\z - z\\ \\z\\ - 2^y^i\\z - z\\ 
= (z- z) T F{z) + fi\\z - z\\{\\z - z\\ - p|| - 2^n) 



(3.9) 



by (3.4). It follows from A = cl that ||x|| i = a/c||x|| , which implies (3.8) by (3.7) and 
(3.9) because \\z — z\\ > 2 ^ > \\z\\ + H= > \\z\\ + 2y/n. Therefore, a zero curve j a of p a is 



< M} by Theorem 



bounded for each a G R 2n with a + A- x F^{a) G B = {z G R 2n 

(2.2), and there exists a 7 a of p a , emanating from (0,a) and reaching a zero z of F^ at 
A = 1 for almost every a G R 2n with a + A~ x F^ x \a) G B. Then, z is a zero of F because 
F^ l \z) = F°(z) = f(z) by Q and Q. □ 



Lemma 3.2. ([4]). For every fi, ipn(a, b) = a + b — ^ (a — b) 2 + 4/i 2 = d if and only if 



A x 
A 2 

p a (A, z) = if and only if 



Let A G (0, 1) and A 



Then, it follows from (3.6) that for a 



a 
a" 



y = f(x) + r, r = fix — (1 — A) [/(a') — a" + pa' — Ai(x — a')], 
v = -jiy + (1 - A)[$ M (a) + pa" - A 2 (y - a")] 
By Lemma 3.2, it follows from (3.10) that p a (A, z) — if and only if 



,1 ; 



2 



> 0, 



-){Vi 



Vi 

2 



p 



(3.10) 
(3.11) 
(3.12) 

(3.13) 



Lemma 3.3. Let a 



a 
a" 



> 



where (3 is the constant given in (|3.5|). Then, we 
have y > for every (X,z) G p _1 (0) with z = (x,y). 

Proof. Let's prove that $ M (a) > . If otherwise, there exists an i such that ip^a^a") = 
< + a "-V( a i -<) 2 + 4// 2 < . Then, a^ + a'/ < ^(a^ - a'/) 2 + 4p 2 and it follows a-a? < p 2 . 
Thus, /3 2 < (3 2 by (3.5) and the condition of Lemma, which is a contradiction. Therefore, 
3V( a ) > • Assume that there exists an index % with yi < . Then, it follows that Vi > 
by gig ) and $ M (a) > 0. But, by gTToJ, 

Ui = ipp(xi,yi) < <p{xi,yi) = 2mm{xi,yi} < 2y l < 0, 

which is a contradiction. Thus, we have y > . □ 



Lemma 3.4. a 



> 



, where (3 is the constant given in (3.5) and (X,z) be 



any limit of a smooth zero curve j a C p 1 (0) . Let 

Pi = {i\xi = oo}, P 2 = {i\yi = oo}, and {i\xi 



-oo}. 



Then, A < 1 implies that i 3 P 2 and fi(x) = — oo for every % G P\ and P 2 = N. And A = 1 
implies that ( i ) Xj < and fi(x) = oo for every i £ P 2 , (ii) fi(x) = oo /or even/ i E N and 
(Hi) yi — and < /or even/ i £ P% . 



^p,(z) and 



Proof. 1) Consider the case of A < 1. If z G Pi, then where v 

p = p(A). Suppose the contrary. Then, there is a finite C > such that Xi — ~ = C , and 
Vi = oo because xi = y + C = oo . Thus, = — oo by (3.12). But, by Lemma 3.3, y% > 
which is a contradiction. Hence, by (3.13), yi — y = 0. Suppose that i E P 2 . Then yi = oo 
and = — oo by (3.12), which implies yi — ^ = oo, contradicting yi — ^ = 0. Therefore, 
i 3 P 2 , i.e. I/, is a finite positive, from which it follows that x,i = oo implies fj = oo, and 
fi(x) = — oo by (3.11). Now, let's prove P 2 = N. Let i G P 2 . Then, we have Vi = —oo by 
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oo by (3.13). Hence, i G N 
oo by (3.11). Because > y by 



(3.12) , and tji — y = oo, which implies x« — y = and x 
and P 2 C iV . Let i £ N . Then, Xj = — oo and fj = 

(3.13) , we have = — oo. Therefore, |/j = oo , i.e. 2 G P 2 by (13.121), which implies N C P 2 . 



Thus, P 2 = N. 



2) Let's consider the case of A = 1. In this case, by ( 3.5[ ) , we have /x = 0, and by ( |3. 13| ) 

> 0. (£<-7f)(f7i 



.if. ; 



Vi 

2 



Vi 

2 



0. 



(3.14) 



First, let's prove the proposition (i). If i G P 2 , then t> j < and ?/j — ^ = oo by (3.12), and 

x 4 -^ = (3.15) 

and fj = and f/j = /j(x) = oo 



by (3.14). Consider the case of fj = . Then, Xi = by (3.15 



by (3.11 ). Consider the case of Vi < . Then, Xi < by 



3.15 ), and < and oo = yi < f%{x) 



by 3.11). Second, let's prove the proposition (ii). If i 6 iV , then u f = — oo by ( 3.14[ ) and 
= oo by ( 3.12[ ), which implies that fi(x) = oo by the proposition (i). Finally, let's prove 
the proposition (iii). If i G Pi , we have Xi — y = oo by (3.12) and 



Vi 







(3.16) 



by (3.14). Suppose that i G P 2 . Then, from the proposition (i), it follows that X{ < , 

by (3.12) because — and A = 1. Thus 
by (3.16), and we have fi(x) + fj = and f j > by (3.11 ). Consider the case of f j = 

> 



contradicting z G Pi . Hence, j 9 ? 2 and f 
. Then |/ 



□ 



fi(x) = 0. Consider the case of > . Then /j(x) + r. 

Theorem 3.2. a=^ a /; ^>^^^, where (3 is the constant given in (|3.5|) and suppose 
that a map f : R n — > R n satisfies /(a') > and 



( X -yf(f( X )-f(y))>Q 



(3.17) 



for some 5 > and for ever y x and y with \\x — y\\ > 5 , and the assumption 1 is satisfied 
for F^ x \z) defined by (3.3). Then, there is a smooth zero curve 7 a C p _1 (0) ; emanating 
from (0, a) and reaching a NCP solution z = (x, y) with y — f(x) at A = 1 . 



Proof. If conditions of the Theorem are satisfied, conditions 
satisfied by Remark 



(iii) of Theorem 2.1 



2.1 



are 



Then, there exists a smooth zero curve 7 a C p _1 (0) of p a (X,z), 
along which the Jacobian matrix Dp a (X,z) has rank n , emanating from (0,a) and not 
intersecting itself by Lemma 2.2 of [4]. Suppose that 7 a is unbounded. Then there is an 
unbounded sequence {(Xk,z k )} C 7 a such that A^ — > A G [0, l],z k — > z and ||z|| = oo as 
k — > oo , where z = (x, y) . Let A < 1 . Then A = is impossible because the zero curve 
7 a can't return the starting point (0, a) and the p a (0, z) has a unique simple zero z = a. Let 
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% 3 F = {i\zi is finite} . First, consider the case of i G P\ . Then, j%{x) = — oo by Lemma 



3.4 Thus, we have 



i X j 



-oo. 



OUiix) - hip!)) 
Second, consider the case of i £ N . Then, x, = — oo and iji 



(3.18) 



oo by Lemma 3.4, and it 



follows from (3.11) that fi(x 



oo . Therefore, we obtain the (3.18) again. Third, consider 

and yi 



the case when Xj is finite. Then jji should be infinite because of i 3 F 
P 2 = N by Lemma 3.4 and Xj = 



oo 



But, 



(]3.18|) for every i 3 F . 

x 



— oo, contradicting the finiteness of Xj . Hence, we have 
'" < 5 for every 5 > 0. Then x — a' because 5 is arbitrary. Hence, x 



Suppose that \\x — a 
is bounded, which implies that y should be unbounded. Since P 2 



N , if yi — oo for some 



i , then Xj = — oo, contradicting xi = a^. Therefore, there is 5 > such that \\x — a'\\ > 5 . 
Then, by (3.17), we have (x — a') T (/(x) — f(a')) > . Thus, taking account of (3.18), there 



is j G F such that 



( x i ~ a 'i)(fi( x ) ~ fi( a ')) 



The following two cases are possible. 



Xj Ctj . 



oo. 



oo, 



(ii) xj < a' p fj(x) 



-oo. 



(3.19) 
(3.20) 



Consider the case (i). Because of j 



(ii). Taking accou nt o: 
Vj > by Lemma 



3.3 



G F 

= u 



is bounded and so is fj by (3.11). Then, we 
+ fj, contradicting j G F. 



x) 



have jjj = oo by (3.19) because jjj = fj(x) + fj, contradicting j G F. Consider the case 

(3.20), we have jjj = — oo similarly to the proof of the case (i), while 
We have a contradiction again. Therefore, we have A = 1 and the 



equality (3.18) for i G P 2 U iV by Lemma 3.4 For i G Pi , we have the equality (3.18) again 



because f(a') > . Hence, it holds (3.18) for every i 3 F . Suppose that ||x — a'\\ < 5 for 
every 5 > 0. Then x = a' , which implies that y is unbounded. Thus, P 2 and then 
Xi < for i G P 2 by Lemma 3.4, contradicting x« = > 0. Therefore, there is 5 > such 
that ||x — a' || > 5 , and we have a contradiction again like the case of the A < 1 . Therefore, 
the zero curve 7 a is bounded and it's accumulation point has form of (1, z) by Lemma 2.3 of 
[4]. Then, x is a solution of (1.1) and y = fix) , where z = (x,y) . The proof is completed. 

□ 



The condition (3.17) for map / : R n — > R n is called generalized monotonicity. 



4 Preliminary numerical experiments 

We can parameterize zero curve 7 a of p a by its arclength s , that is, there exist continuously 
differentiable functions x(s) and A(s) such that 



g a (x(s), A(s)) = 0, x(0) = x°, A(0) = 0. 



(4.1) 
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By differentiating the first equation of (4.1), we obtain the following result: the homo- 



topy path is determined by the following initial value problem of the ordinary differential 
equations. 

D ea (x(s),\(s))( X A = 0, 

x(0) = x°,A(0) = 0, ( 4 - 2 ) 
|| (±(s),A(s)) || =1,A(0) >0,, 

where Dg a (x(s), A(s)) = (^f, ^f) ■ Then the x component of the solution (x(s*), X(s*)) 



of (4.2), which satisfies A(s*) = 1, is the solution which we need to find. We carried out 
numerical experiments for nonlinear equations to show the performance of our NFPH method 
by PC Pentium IV, 3.19GHz, 1.00GB of RAM. 

Our homotopy zero-curve trace was made by using ODE (ordinary differential equation) 
toolbox of MATLAB. In our experiment, we took A = oel, a > and made use of 'ode45' or 



'odel5s', ODE solvers of MATLAB, to solve (4.2). 



Let Sf be the upper bound used in the ODE part for arc lengths and C n be the 
number of intermediate checks Si,S 2 ,... between and Sf. Then interval [0,5^] is di- 
vided in (C n + 1) intervals of equal length. The zero-curve trace is finished as soon as 
one candidate solution has been found in any of the intervals defined by the checkpoints 
[0, Si, S%, Sf]. If path-following is successful, we obtain the more refined solution 'nsol' 
using 'fsolve', nonlinear equation solver of MATLAB, with initial guess 'hsol' obtained by 
homotopy method. In the following tables, a is the starting point x°, sol is a solution of the 
given problem, N c is the number of checked intervals, and fhom = F(hsol)/ (1 + ||/iso/||), 
fnew = F{nsol)/ (1 + ||nso/||). 

Example 1. 



F(x) = 2x - 4 + sin(2nx),x E [-100, 100], a = 0, sol = 2. 

Example 2. 

F(x,q) = f rt^Zg) ) E [-100,100] x [-100, 100], a = [0, 0] 

Example 3. 

/ x + 0.5g + 0.3^-5 \ 
F(x,q,z)=l 0.6x + q + O.lz -7 

\ 0.2x + OAq + z - A J 

{x,q,z) E [-100,100] x [-100,100] x [-100, 100], a = [0,0,0] 

Example4. 

F(x) = arctan{100x)/-K + sm(5x/(x2 + 0.2))/2 + 0.1x,x G [-2, 2], a = 0.2, sol = 0. 
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The results of numerical experiments are shown in the following tables. 



Tablel(Examplel, Sf = 2.5, C n = 70) 



method 


N c 


hsol 


nsol 


time(s) 


a = 0.001 


14 


2.0005 


2.00000 


1.0809 


a = 50 


2 


1.9783 


2.00000 


0.3327 


FPH 


20 


2.0008 


2.00000 


1.5547 


NH 


14 


2.0000 


2.00000 


0.9153 


Table2(Examplel, Sf = 


z 5, C n = 


70) 


method 


N c 


hsol 


nsol 


time(s) 


a = 0.001 


7 


2.0005 


2.00000 


0.5911 


a = 50 


1 


2.0032 


2.00000 


0.1862 


FPH 


10 


2.0008 


2.00000 


0.7333 


NH 


7 


2.0000 


2.00000 


0.5098 



Table?,{Example2, S f = 20, C n = 70) 



method 




hsol 


nsol 


fhom 


fnew 


time(s) 


a = 0.001 


3 


-7.4234e-001 
-6.7601e-001 


-7.3908e-001 
-6.7361e-001 


4.0222e-003 
-2.7211e-006 


3.5192e-011 
3.6061e-012 


0.5689 


a = 50 


1 


-7.2803e-001 
-7.3465e-001 


-7.3908e-001 
-6.7361e-001 


3.4288e-003 
3.4037e-003 


1.6008e-012 
1.7203e-013 


1.4577 


FPH 


6 


-7.3885e-001 
-7.3465e-001 


-7.3908e-001 
-6.7361e-001 


1.2271e-003 
1.1222e-003 


2.3637e-012 
5.1031e-013 


1.2431 


NH 


3 


-7.3571e-001 
-6.7112e-001 


-7.3908e-001 
-6.7361e-001 


-4.1720e-003 
-4.5865e-003 


4.1332e-011 
4.2909e-012 


0.3311 



Table4(Example2, Sf = 5, C n = 70) 



method 




hsol 


nsol 


fhom 


fnew 


time(s) 


a = 0.001 


10 


-0.71376658 
-0.65472611 


-0.73908522 
-0.67361213 


-3.14294e -02 
2.059234e-05 


1.3983767e- 07 
1.4943515e-08 


1.4112 


a = 50 


1 


-0.72803942 
-0.73465009 


-0.73908513 
-0.67361202 


3.428826e-02 
3.403779e-02 


1.6008305e-012 
1.7202905e-013 


1.4576 


FPH 


21 


-0.73915876 
-0.67293244 


-0.739085321 
-0.673612169 


-4.03203e-04 
-3.67077e-04 


2.337165e-07 
9.0814323e-010 


3.5987 


NH 


10 


-0.7142857 
-0.6550778 


-0.739085223 
-0.673612122 


-3.080903e-02 
-4.510372e-016 


1.288977e-07 
1.378482e-08 


0.8984 
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Table5(Example3, S f = 30, C n = 50) 



method 


N c 


hsol 


nsol 


fhom 


fnew 


time(s) 






1.67897430 


1.67155427 


3.0576e-003 


4.3102e-012 




a = 0.001 


3 


5.89119023 


5.86510265 


4.2828e-003 


9.4505e-012 


0.9735 






1.32550592 


1.31964807 


2.4461e-003 


1.5067e-011 








1.66456368 


1.67155425 


2.8405e-003 


-1.5842e-011 




a = 50 


1 


5.92393350 


5.86510263 


7.4144e-003 


-3.3679e-011 


0.9152 






1.31394188 


1.31964809 


2.2530e-003 


-1.3984e-012 








1.6731148 


1.6715542 


1.7504e-003 


-1.3416e-011 




FPH 


2 


5.8574770 


5.8651026 


2.4505e-003 


2.0760e-012 


0.6623 






1.3209060 


1.3196482 


1.4003e-003 


1.0126e-011 








1.6758000 


1.6715542 


1.7504e-003 


-1.3416e-011 




NH 


3 


5.8800000 


5.8651026 


2.4505e-003 


2.0760e-012 


0.3325 






1.3230000 


1.3196482 


1.4003e-003 


1.0126e-011 





Table6(Example3, S f = 2,C n = 50) 



method 


N c 


hsol 


nsol 


fhom 


fnew 


time(s) 


a = 0.001 


33 


1.50807356 
5.29045112 
1.19058642 


1.67155426 
5.86510262 
1.31964809 


-0.07385090 
-0.10343821 
-0.05908049 


1.680643e-011 
-1.5938702e-009 
-6.454418e-010 


6.7216 


a = 50 


1 


1.69024906 
5.79197983 
1.33466574 


1.671554252 
5.865102639 
1.319648094 


-0.001861054 
-0.008413507 
-0.001461469 


6.8498442e-012 
2.679456e-012 
5.0761474e-011 


0.7997 


FPH 


26 


1.7520000 
5.3961542 
1.3848495 


1.67155424 
5.86510263 
1.31964809 


-0.01965900 
-0.06054967 
-0.01553925 


-4.948133e-010 
-4.022569e-010 
-3.449198e-010 


4.6629 


NH 


33 


1.5047999 
5.2800000 
1.1880000 


1.67155425 
5.86510264 
1.31964809 


-0.075378066 
-0.105529292 
-0.060302452 


6.0140810e-010 
4.7703930e-010 
6.6359281e-010 


1.6579 



T able! [Example^ Sf = 5, C n = 70) 



method 


N c 


hsol 


nsol 


time(s) 


a = 0.001 


37 


-1.8652e-005 


3.1351e-011 


4.8317 


a = 1 


18 


1.9012e-004 


7.6041e-020 


2.3264 


a = 75 


1 


-2.3974e-004 


5.3332-022 


0.3118 


FPH 


7 


-6.0612e-006 


1.0732e-012 


0.8831 


NH 


37 


5.0773e-005 


2.9810e-022 


3.3849 



As shown in Table 1 and Table 2, when a = 50, our NFPH is much better than FPH and 
NH for Exampl 1. For Example 2, our NFPH with a = 50 is better than others in the 
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accuracy of obtained solution (Table 3 and Table 4). Our NFPH (C/ = 2, a = 50) is 
much better than others in both of time and accuracy(Table 6) for Example 3. Especially, 
the NFPH with a = 75 shows good performance for Example 4 (Table 7). From the 
above numerical results, we can see that choice of proper A can remarkably improve the 
performance of the NFPH. 

5 Conclusion 

This paper describes a probability-one homotopy algorithm for solving nonlinear systems 
of equations and complementarity problems. They are attractive because they are able to 
solve a qualitatively different class of problems than methods relying on merit functions. 
This claim is justified both theoretically and computationally. While the common homotopy 
used to solve nonlinear system F(x) = is FPH defined by p(a, A, x) = XF(x) + (l — X)(x — a) : 
in this paper we considered probability-one global convergence of the algorithm based on 
NFPH defined by 

p(a, A, x) = \F(x) + (1 - \)(F(x) - F(a) + A(x - a)) 

and extended the results to NCP with generalized monotonicity. The preliminary numerical 
experiments for some difficult nonlinear equations showed the robustness and fast conver- 
gence of the NFPH method. We expect the NFPH method would have better performance 
than NH(Newton Homotopy) or FPH method because NFPH combines the advantages of 
both NH and FPH. 
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